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' Abstract 

oo : 

T^ij- , Two cosmological models with non-phantom matter having the same expansion 

CN I of the universe as phantom cosmologies are constructed. The first model is charac- 

terized by the evolving gravitational "constant" G and a dark energy component 
(3 ■ with a non-conserved energy-momentum tensor. The second model includes two 

I interacting components, the dark energy component and the matter component. 

Qh' Closed form solutions are obtained for the constant values of model parameters 

Q ■ and constraints on the parameters of each model from cosmological observations 

are outlined. For both models it is explicitly shown how the components of each 
^ I model produce the expansion of the universe characteristic of phantom cosmologies, 

despite the absence of phantom energy. These findings stress the interpretation 
of phantom energy as an effective description of the more complex dynamics of 
^ ' non-phantom matter. 



1 Introduction 

Complementary cosmological observations of supernovae of the type la (SNIa) pQ, 
cosmic microwave background radiation (CMBR) 2^, large-scale structure (LSS) 
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and other cosmic phenomena have firmly estabhshed the picture of the accelerated 
expansion of the contemporary universe _4 . The present accelerating phase of the 
expansion of the universe and its onset at a relatively low redshift {z 1) represent 
one of the most intriguing and most studied problems in modern cosmology. The 
majority of theoretical explanations of this phenomenon invoke a new component 
of the universe named dark energy ^ with the equation of state (EOS) 

Pd = wpd , (1) 

where pd and pd represent dark energy density and pressure, respectively. The 
most studied theoretical candidates for the role of dark energy are the cosmolog- 
ical constant (CC) {w = —1) [HI 13 IE), and its dynamical variants such as the 
renormalization group running CC [HlEl^B! quintessence {w > —1) ^j, tachyon 
models ^31 (w > —1) and the Chaplygin gas {w > —1) 

Recent analyses of cosmological observations ^1 El El El E] allow, and even 
favour, a sort of dark energy with a supernegative EOS, i.e. w < —1. This un- 
orthodox type of dark energy, first introduced in [2U| , was named phantom energy. 
Many models of phantom energy appeared soon |2^, addressing both its funda- 
mental implications and cosmological consequences. One of the most interesting 
features of phantom energy is certainly the possibility of the divergence of the scale 
factor of the universe in finite time. The expansion of the universe in such a model 
of phantom energy leads to the unbounding of all bound structures, a phenomenon 
also vividly referred to as "Big Rip" 22j. 

Although phantom energy represents a phenomenologically appealing possibil- 
ity, the violation of the dominant energy condition (DEC), inherent in phantom 
energy models, leads to problems at the microscopic level. For example, it is pos- 
sible to describe phantom energy in terms of the effective scalar field theory with 
negative kinetic terms, valid up to some cut-ofF scale. In such a formulation the 
vacuum of the theory is no longer stable, i.e. phantom energy decays. Such theories 
can still be cosmologically viable if the lifetime of phantom energy surpasses the 
age of the universe. This requirement puts stringent constraints on the parameters 
of the effective scalar filed theory, above all on its cut-off scale [23^ i24j . There still 
remains a question whether some other viable microscopical formulation exists. 

In such a conflict between favour from the observational side and disfavour 
from the theoretical side, phantom energy models face an interesting alternative: 
possibility that matter which has no phantom characteristics (e.g. satisfies DEC) 
produces observational effects attributed to phantom energy. In this paper we con- 
sider two realizations of this possibility. The first realization given in section [21 is 
a model reminiscent of generalized phantom energy j25j . in which we consider a 
sort of cosmology with a time-dependent gravitational "constant" G and a dark 
energy component with a non-conserved energy-momentum tensor. The second 
realization is based on the dynamics of two interacting cosmological components 
and is displayed in section |31 

^ There are alternative explanations rooted in brane- world models which do not require dark energy 

u 
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2 A model with an evolving G 



We consider a cosmological model with two components. The first component, 
which we caU the matter component, has the equation of state 

Pm = 'l{a)pm , (2) 

where 7(a) > 0, and p„i and pm denote the energy density and pressure of the 
first component, respectively. We assume that the energy-momentum tensor of 
this component is conserved, T^^^ = 0, which leads to a well-known relation for 
the scaling of pm with the scale factor a ^: 

Pm = Pmfl e -^"O . (3) 

The second component, which we call the dark energy component, satisfies DEC 
and has the equation of state 

Pd = v(.a)pd , via) > -1 ■ (4) 

Here pd stands for the energy density and pd denotes the pressure of the dark 
energy component. We assume that the energy-momentum tensor of the dark 
energy component is not conserved, i.e. T^^ 7^ 0. Therefore, the parameter of EOS 
(jlj) does not determine the scaling of pd with a. One possible way of harmonizing 
non-conservation of Tj^'^ with the general covariance of the Einstein equation is 
the promotion of the gravitational constant G into a time-dependent function G{t) 
(see reference [251 details) ^. G{t) satisfies a generalized conservation condition 

(G(t)r^n;. = 0, (5) 

where T^'^ = + T^'^ . It is important to stress that the procedure explained 
above does not represent some trivial multiplication of the constant G by some 
function of time f{t) and multiplication of the total energy-momentum tensor T^^^ 
by f{t)~^ since the energy-momentum tensor of the matter component is conserved. 
The relation © can be expressed as 

d{Gpd) + PmdG + ZGpdil + r?(a))— = . (6) 

a 

This equation determines the dynamics of G in terms of energy densities and 
parameters of EOS of the components of the universe. At this place, it is important 
to notice that in the Friedmann equation 

d\2 /c _ Svr 
a) a? 3 



7) +-:2 = —'^iP^^ + Pd)^ (7) 



the evolutions of both G and energy densities pm and pd determine the kinematics 
of the universe, i.e. the function a{t). The aim of this section is to investigate 



^The subscript denotes the present epoch throughout the paper. 

^Many models consider the time-dependent G, such as the renormalization group running of G 
ini El Uni J the time-dependence of G originating from extra dimensions [221 or the effective G in 
scalar-tensor theories |28| . 
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the possibility that the product G{pd + Pm) has a component which grows with 
the scale factor (i.e. its effective parameter of EOS is smaller than —1), while the 
components of the universe have non-phantom nature (they satisfy DEC). 
To this end, we introduce an assumption of the following scaling behaviour: 



Gpd — Gopdfi 



(8) 



where w{a) < —1. This assumption clearly introduces a source into the Friedmann 
equation ((7| which is identical with the source originating from phantom energy 
with the parameter of EOS w{a) in a model with constant G. 
The evolution equation for G becomes 

-3{l+w{a))-l 



dG 



^GfiPdfl 



dw{a) + {w{a) 



"^In 


r a \ 






.00 


yaoJ 



V{a))d(- 



(9) 



Generally, it is not possible to solve this equation in closed form, so further in this 
section we consider a simplified model with 7(a) = 7 = const, ri{a) = r] = const 
and w{a) = w = const, to gain deeper insight via an analytical solution which one 
can obtain in this case. The function G can now be expressed in terms of a as 

—3(111—7) 

G = Go I 1 



Pdfi rj-w 



Pmfl -y -w 



1 



(10) 



Once we have the expression for G, we can give the expression for the other source 
term in the Friedmann equation (the first is given by (jHI): 

-3(1+7) n — w / " ^ ~3(l+ii') 



Gpm — Go Pm.,0 + Pdfi 



V 



W 



GoPdfl 



rj 



w 



(11) 



7 — w/Vao/ 7 
while the total source term (the right-hand side of the Friedmann equation) be- 



comes 



G{pm + Pd) = Go Pm,0 + Pdfi 



rj — w 

) 

J — w 



-3(1+7) 



+ GoPd, 



1-V 



7- 



w 



(12) 

Closer inspection of equation (fT^ shows that in our cosmological model the uni- 
verse evolves as if the quantity G were constant and we had one phantom compo- 
nent with the parameter of EOS w and one non-phantom component with the pa- 
rameter of EOS 7, both components having conserved energy-momentum tensors. 
Therefore, our cosmological model mimics the behaviour of the model with a non- 
phantom matter component with the present energy density pmfi = Pmfi + Pdfi 
and a phantom component with the present energy density pdfi = Pdfi^_w 
acceleration of the universe in our model is given by 

-3(1+7) 

'1 + 37) ( Pmfi + 



r]—w 
7— ui 

^ The 



a 

a 
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+ 1 + 3r/ - (1 + 37) 
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Let us finally consider possible constraints on the parameters of the model. The 
assumption of constancy of all parameters of EOS is probably an oversimplification. 
However, for modest variations of some of the parameters with the scale factor a, 
one expects that the model given above represents a good approximation. It is 
certainly conceivable that for the more realistic parameters w{a), rj(a) and 7(a) one 
can, at the level of the Priedmann equation, obtain in our model a sort of dynamics 
which is identical with the dynamics in the more general model with constant 
G and conserved energy-momentum tensors of the phantom and non-phantom 
components. Given many observational constraints on the past evolution of G j29j . 
one would expect quite stringent constraints on the difference rj — w, which should 
be small. This means that this model could explain the cosmological expansion 
with w not much more negative than —1. However, in the general case of our model, 
a sort of cosmology where w would differ from — 1 more substantially is certainly not 
excluded. The expression for the acceleration of the expansion of the universe 
provides another constraint on the parameters of the model. Namely, in order 
to have a transition from deceleration to acceleration at low redshift (z ~ 1), the 
coefficient l+3ry— (1+87)^5^ must be negative. In the case that the difference rj—w 
is small, this requirement reduces to the standard one, r] < —1/3. The requirement 
oi rj — w being small also favours the possibility rj = —1, which is equivalent 
to the time-dependent cosmological constant. Models with the time-dependent 
cosmological constant and G, studied in represent a specially interesting 

case. Models with a growing cosmological constant (and a time-dependent G) 
exhibit a very peculiar fate of the universe, leading to the unbounding of all 
gravitationally bound systems, while leaving non-gravitationally bound systems 
unaffected, the so-called "partial rip" scenario. 



3 A model with two interacting components 

In this section we consider a model with two interacting, non-phantom components. 
The first component, the dark energy component, is described by the equation of 
state 

Pd = il{a)pd, r/>-l, (14) 

while the other component, the matter component, is determined by the following 
equation of state: 

Pm, = 7{a)Pm , 7 > . (15) 

In this model, the gravitational constant G has no space-time variation. The 
interaction of the components is included in the model in the following way. We 
assume that the energy- momentum tensors of the two separate components are not 
conserved, but the total energy- momentum tensor T^'' = + T^'^ is conserved. 
In this way, there exists an exchange of energy and momentum between the two 
components. The requirement of the conservation of the total energy-momentum 
tensor can be expressed as 

dpm + 3pm{l + 7(a))— = -dpd - 3pd(l + r]{a))— . (16) 
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What remains to be determined is the specification of the interaction (energy- 
momentum exchange) between the components. The aim of this model is to demon- 
strate that this set-up can mimic the expansion of the universe characteristic of 
phantom cosmologies. Therefore we assume the following evolution law for the 
dark energy component: 

^ ^ ~3{l+w{a)) 



Pd = Pd,o — , (17) 

where w{a) < —1, i.e. the dark energy component has the evolution law char- 
acteristic of phantom energy. The non-phantom dark energy component has the 
scaling with a characteristic of phantom energy owing to the interaction with the 
matter component. Equation H16() then determines the evolution law for the en- 
ergy density of the matter component. For general values of the parameters of 
EOS it is not always possible to obtain the solutions in closed form. Therefore, in 
the remainder of this section we assume that these parameters are constant, i.e. 
7(a) = 7 = const, r]{a) = rj = const and w(a) = w = const. This particular choice 
will allow us to gain insight via closed form solutions. The energy density of the 
matter component then becomes 

r/-u;\ / a\-3{i+7) ry - w; / a \ ^^(1+"') , , 

Pmfl + Pd,o — - Pd,o — • (18) 

J -wj \ao/ 7 - w \aoJ 

The total energy density, p = Pm + Pd, which appears on the right-hand side of 
equation 0, then has the form 

/ T/-u;\ / a \-3{i+7) ^_ri / a \ -^(1+"') , , 

P={Pm,o + Pd,o — Pd,o{ — ] ■ (19) 

V -wj \aoJ -w \ao/ 



The acceleration of the expansion of the universe is given by the expression 

- = (1 + 37) Pmfl + Pdfl- I , 

a 6 \ 7 — u)/\ao 

3(1+M«)' 



r]-w \ / a \ -3(1+7) 



/ , ^n — w\ (a 
+ l + 3r?- 1 + 37-^ ]pdfl[- 



(20) 



As for the model displayed in section 121 the right-hand side of the Priedmann equa- 
tion is the same as in a model with constant G and two non-interacting compo- 
nents: the first being the non-phantom component with the present energy density 
Pmfl = Pmfl + Pdfl^p^ and the parameter of EOS 7, while the second being phan- 
tom energy with the present energy density pdfl = '^E^Pdfl and the parameter of 
EOS w. Equation ()18|) shows the effects of the interaction with the dark energy 
component on pm as an additional term growing as a~^^^~^^\ The requirement 
that the scaling law of the matter component should not differ too much from 
the scaling law dictated by its EOS (~ a^^^^^"'^), i.e. that the interaction is not 
too strong, leads to the condition that the difference r] — w should be small. In 
the model with a more general variation of some of the parameters 7, ij or w, it 
is conceivable that this constraint would be milder. Again, as in section 121 two 
non-phantom components mimic phantom cosmology. The model can successfully 
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describe the transition from the decelerating to the accelerating regime of the ex- 
pansion of the universe if the coefficient 1 + 3ry — (1 + 87)^5^ is negative. When rj 
is close to w, the afore-mentioned requirement reduces to the condition rj < —1/3. 
One especially interesting variant of the model is the case r] = —1. In this case, the 
evolving cosmological constant in interaction with the matter component mimics 
the expansion of phantom cosmology. 

One way of elaborating the model given in this section would certainly be its 
formulation in terms of classical fields. We can model the system of two interacting 
components as a system of two minimally coupled interacting scalar fields in a 
cosmological setting. For the Lagrangian of the interacting system we then take a 
general form (we consider only time-dependent scalar fields) 

^ = Y + Y-m^), (21) 

where cf) and V denote scalar fields. Given that the total energy density \s p = 
^ + ^ + V{(t), ip) and the total pressure is p = C, one obtains the following two 
constraints on the dynamics of the scalar fields: 

</,2+^2 ^ (1 + + (1 + ^)^^ ^ 

2V{^,ij) = {l-r,)pd + {l-l)Pm. (22) 

From (fTTj) and (fTH|) we have obtained pd and pm , respectively, as functions of the 
scale factor a. On the other hand, from ((71) we can determine the function a{t). 
This makes the right-hand sides of equations (|22|) known functions of time. All 
pairs of the functions (p and ip (with a nontrivial potential V{(j),ip)) that satisfy 
equations (|22)) can produce the evolution of the universe as described in the model 
of this section. This class of solutions certainly does not exclude more sophisticated 
(and realistic) field (or microscopic) models. 

4 Conclusions 

The two models, described in sections El and |S1 have been constructed to demon- 
strate that cosmologies without phantom energy can lead to an expansion of the 
universe usually attributed to phantom energy. The first model is characterized 
by an evolving gravitational "constant" G and a dark energy component with a 
non-conserved energy-momentum tensor. The second model is based on two in- 
teracting components. Both models yield results for the cosmological evolution 
of their components which are testable against the results of various cosmological 
observations. Calculations in this paper have been made with a specific choice of 
parameters (e.g. constant parameters of EOS) which ensures closed form solutions. 
These solutions facilitate the interpretation of the physical meaning of the obtained 
results, but the scope of the models described in this paper certainly does not end 
here. Models with variable (e.g. dependent on a) parameters 7, r] and w offer much 
more possibilities (especially in terms of satisfying numerous constraints from the 
past evolution of the universe) and merit further investigation. The possibility of 
mimicking phantom cosmology by non-phantom one, certainly does not rule out 
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an appealing and provocative idea of phantom energy. However, it puts a greater 
ponder on the nature of phantom energy as an effective description of the more 
complex dynamics of non-phantom matter. 
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